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Abstract: Pearson's R is the most common correlation statistic, used mainly in parametric set- 
tings. Most common among nonparamctric correlation statistics are Spearman's S and Kendall's T. 
We show that for bivariate normal i.i.d. samples the pairwise asymptotic relative efficiency (ARE) be- 
tween these three statistics depends monotonically on the population correlation coefficient p. Namely, 
ARE H>T (|p|) increases in |p| £ [0, 1) from 1.096 ... to 1.209 . . . , ARE fljS (|p|) increases from 1.096 . . . 
to 1.439... and AREt,s(|p|) increases from 1 to 1.190.... This monotonicity is a corollary to a 
stronger result, which asserts that q R ,T-.a(p) = (ARE^x(p) — AREjj t(o) — ARE^ T (a)(p— a))/ (p— a) 2 
is increasing in p £ (0, 1), where a £ {0, 1}, and similarly for the two functions q R g. a and qx s-a with 
ARE/j s and AREt,s m place of ARE/j t- Another immediate corollary is the existence of quadratic 
(in p) polynomials L R}T . a (p) and U R „T;a(p) such that L R,T;a(p) =S ARE J?>T (p) < U R ^ T . a {p) for 
p £ (0, 1) and a £ {0, 1}, and similar quadratic bounds are given for ARE/^5 and AREy ^. The 
proofs rely on the use of PHospital-typc rules for monotonicity patterns. 
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1. Introduction 

Pearson's R, Spearman's S and Kendall's T are the three most commonly used correlation statistics, the 
latter two especially in nonparametric studies. When the population distribution is bivariate normal, the 
question of independence between the two random variables (r.v.'s) reduces to deciding if the population 
correlation p is 0. In the case of testing H : p = 0, it is known that the Pitman asymptotic relative efficiency 
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(ARE) of R to S is ^ [8] and that of T to S is 1 [13] (and hence the ARE of R to T is ^ as well). While 
perhaps less common in practice, one could also use any three of these statistics to test hypotheses of the 
form H : p = p a (against alternatives p > p , p < p , or p =/= p ) f° r arbitrary p G (—1,1). In [2], values 
of the AREs,r(po) (the ARE of S to R for the null hypothesis p — po) are tabulated for several values of 
po € [0, 1); several values of KREt,r{po) are given in [14] as well. 

In this paper, we show that ARE^-rO^ol) is strictly increasing in |po| € [0, 1) from 1.096 ... to 1.209 . . . , 

AREfl ; s(|p |) increases from 1.096... to 1.439..., and AREt,s(|po|) increases from 1 to 1.190 Thus, 

all these ARE's stay rather close to 1 for all values of p € (—1,1). Additionally, several upper and lower 
quadratic bounds are shown to take place for each of ARE/^t, ARE^s and AREt.s- All of these results are 
immediate corollaries to a stronger result, stated in this paper as Theorem 2.1. 

For testing H n : 9 = 9 in the framework of a given statistical model (against any of the alternative 
hypotheses 9 ^ 6$, 9 > 9q, or 9 < 9q), under certain general conditions there exists an easily applicable 
formula for computing the ARE between two (sequences of) real-valued test statistics T\ — (Ti, n ) ne ^ and 
T2 = (T 2 , n )n£N- The main condition (see e.g. [16, 7, 11]) is that the distribution function (d.f.) of cither 
properly normalized test statistic converges to the standard normal d.f. <f> uniformly in a certain sense as 
the sample size n tends to 00. Particularly, if there exist continuous real- valued functions ptTj and on 
the parameter space 6 such that 

supsup pj ^zmp. ^ z \. m _ , 

where V is some neighborhood of 6*o chosen such that pt, is continuously differentiable and ctt, > on V 
for j = 1,2, then the ARE of T\ to T 2 may be expressed by the formula 

<4 2 (<?o) A^o) 2 
ARE(6» ) := ARE Ti ,t 2 (0o) = , a \ 1 1 m , 2 , 

assuming that p,' T .(6o) > 0. The functions and cftJ^ti may be called the asymptotic mean and standard 
deviation, respectively, of the sequence Tj. 

Berry-Esseen bounds provide a nice way to verify the condition (1.1). Such bounds for the Kendall and 
Spearman statistics, which are instances of so-called U- and ^-statistics, are essentially well known; see 
e.g. [12]; in fact, we are using here a result by Chen and Shao [3] and a convenient representation of any 
^-statistic as a [/-statistic [6]. As for a Berry-Esseen bound for the Pearson correlation statistic, we are 
using an apparently previously unknown result in [20]. 

According to the formula (1.2), the ARE between two test statistics can be expressed in terms of the 
asymptotic means and variances of the two statistics. In turn, the asymptotic variance of either T or S in 
the bivariate normal model can be expressed using Schlafli's formula [22] for the volume of the spherical 
tetrahedron in R 4 . Such formulas have been of significant interest to a number of authors; see e.g. the 
recent papers [10] and [15]. We remark also that Plackett [21] obtained a result more general than Schlafli's. 
Actually, here we are using formulas by David and Mallows [5] which are based on [21]. 

To prove the main result, we use l'Hospital-type rules for the monotonicity pattern of a function r = £ 

on some interval (a, b). Knowledge of the monotonicity of ^7 on (a, b), along with the sign of gg' on (a, 6), 
allows one to obtain the monotonicity pattern of r; see Pinelis [19] and the bibiliography there for several 
variants of these rules and applications to various problems. For convenient reference these rules are stated 
as Theorems A, B, and C in Section 3.2. 



(1.1) 



(1.2) 
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2. Monotonicity properties of the ARE in the bivariate normal model 

Let (V n ) =: ((X n ,Y n )) be a sequence of independent, identically distributed (i.i.d.) nondcgcnerate bivariate 
normal r.v.'s with 



pa X (Ty o Y 



EVi =: (fix , Vy) and Cov(Vi) 
note that p € (—1,1) is the correlation coefficient between X\ and Yi. Let 

HU{Xi-X){Yi-Y) 



(2.1) R := i?„ := 



where (X, Y) :— - Y^i=i Vi> R i s commonly called Pearson's product-moment correlation coefficient, and it 
is the maximum-likelihood estimator of p. Spearman's rank correlation is 

(2.2) S := S n := £ KW) - 3(n± T, 

n 6 — n ' n — 1 

i=l 

where r(Xj) := X,j=i '{^j ^ and : = X,&=i '{^fc *s ^i} arc the ranks (and l{-} denotes the indicator 
function). Note that S is simply the product-moment correlation of the sample of ranks (r(Xi),r(Yi)), . . . , 
(r(X n ),r(Y n )). Let next 

Jij := \{Xj < Xi} \{Yj < Y^, 

and let 

(2-3) T:=T n :=-L £ 

denote Kendall's correlation statistic, where fiT(Vi,Vj) := 2(Jjj + Jji) — 1, so that almost surely (a.s.) 
hr{Vi,Vj) — ±1 depending on whether the pair (Vi,Vj) is concordant or discordant; also, E hriVi, Vj) = 
if Vi and Vj are independent. 

Consider the hypothesis test H Q : p = p a against the alternative H\ : p ^ p (or again, either of the two 
one-sided alternatives), where po € (—1, 1). We shall show that each of R, S, and T satisfies the condition 
(1.1), so that (1.2) may be used to express the ARE between any two of these statistics. Further, it is easy 
to see, and also will be clear from what follows, that u\, and a\ are all even functions of p, and also 
p,R, ps, an d Pt are odd functions, so that the ARE of any pair of these statistics is even. See Figure 1 
for a plot of these three functions, and note it suggests each of the pairwise ARE's is strictly increasing on 
(0,1). Further, the shapes of these plots suggest the functions may be well- approximated by a quadratic 
polynomial. Indeed, the monotonicity of the ARE and a quadratic approximation shall be immediate results 
of the following: 

Theorem 2.1. For the test of the null hypothesis p = p against any of the three alternative hypotheses: 
P 7^ Po, P > Po, or p < p , let 

ARE Tl , T2 (p )-ARE TliT2 (a)-ARE^ iiT2 (a)(po-a) 

9a(Po) := qT u T 2 ;a{P0) : = 7 ^ ' 

(Po - a) 

for po £ [0, a) U (a, 1) and a G [0, 1], where T; t is one of the statistics R, S, or T; here and in what follows, 
ARET!,T 2 ( a ) and ARE^ T2 (a) are understood to mean AREti,t 2 (1 — ) and ARE' Tl T2 (l—) , respectively, when 
a = 1. Then 
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(RTO) q Rt T-,o is increasing from 0.0966... to 0.1125... 

(RT1) q R>T -i is increasing from 0.1510... to 0.2247... 

(TSO) qT,s-.a is increasing from 0.0984- ■■ t° 0.1904- ■■ 

(TS1) 9t,S;1 is increasing from 0.5516... to 1.8200... 

(RSO) q Rj s;0 is increasing from 0.2045... to 0.3428... 

(RSI) q RtS -i is increasing from 0.8682. . . to 2.6639. . . . 

The term "increasing" will mean for us "strictly increasing," and similarly "decreasing" will mean "strictly 
decreasing." Exact expressions for the endpoint values 0.0966. . . , 0.1125. . . , ... are given at the end of the 
respective sections RTO, RT1, ... in the appendices. All proofs arc deferred to Section 3. 
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Fig 1. Plots of ARE r,t(po) (solid), ARE T ,s(po) (dashed) and ARE r,s(po) (dotted) 



Corollary 2.2. For the test of the null hypothesis p = p a against any of the three alternative hypotheses: 
p ^ po, p > p Q , or p < p Q , one has 

(RT) AREjj )T (|po|) is increasing in \p \ G (0, 1) from ^ = 1.0966. . . to = 1.2091 ...; 

(TS) ARE T ,s(|p |) is increasing in |p | € (0, 1) from 1 to ^"jf ~ 15) = 1.1904 

(RS) ARE^sOol) is increasing in \p \ e (0, 1) from ^ = 1.0966. . . to 3(11 ^ 15) = 1.4395. . . . 

This corollary justifies the conjecture that the pairwise ARE's are increasing on (0, 1), which one would 
make from observing Figure 1. 

Corollary 2.3. Let 

L a (x) := L TuT2 . a (x) := ARE TuT2 (a) + ARE' TuT2 (a)(\x\ - a) + g Tl>T2 . a (0+)(|z| - a) 2 , 
U a (x) := U Tl ,T 2 ;a(x) := ARE Tu T 2 (a) + ARE' Ti T2 (a)(\x\-a) + q Tl . T2 , a (l-)(\x\^a) 2 , 
L(x) := L Tl! T 2 (x) := L Ti ,t 2 ;o(x) V L Tlt T 2 -,i{x), 
and U(x) := U Tu t 2 (x) := U Tu t 2 , (x) A U Tii t 2 -i{x) 

for(T u T 2 ) S {(R,T),(T,S),(R,S)}, xG (-1,1), and a G {0,1}. Then for all p G (-1,0) U (0,1) 

La(po) ^ L{po) < ARE(po) < U(po) < U a ( Po ). 
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These piecewise quadratic bounds are illustrated in Figure 2. Note that L and Uq give good quadratic 
approximations to the ARE near the origin, while L\ and U\ are better approximations when p is near ±1. 




(b) Rescalcd plots of /ti,t 2 f° r / = ARE (solid), / = L (dotted) and / = U (dashed) 
Fig 2. Illustration of piecewise quadratic bounds of Corollary 2.3 



Remark 2.4. Numerical approximations to the various bounds in Corollary 2.3 are given below. 



Lr,t-,o{x 
Ur,tA x 

Lt,S;()( X 
Ut,S;0(x 

Ur,S;()( x 



1.0966 + 0.0966a; 2 ; 
1.0966 + 0.1126a; 2 ; 
1 + 0.0984a; 2 ; 
1 + 0.1905a; 2 ; 
1.0966 + 0.2046a; 2 ; 
1.0966 + 0.3429a; 2 ; 



Lr,t-,i{x 
Ur,t-,i{x 

Lt,S;i( x 
UT,S;l(x 
Lr,S;i( x 
UR,S;l{x 



1.0966 - 0.0384|x| + 0.1510a; 2 ; 
1.1704 - 0.1860|a;| + 0.2248a; 2 ; 
1 -0.3612|a;| + 0.5516a; 2 ; 
2.2684 - 2.8980|x| + 1.8200a; 2 ; 
1.0966 - 0.5254|x| + 0.8683a; 2 ; 
2.8924 - 4.1169|d + 2.6640a; 2 . 



Raymond Molzon and Iosif Pinelis/Monotonicity properties of the ARE 



6 



Further, one has 

Lr,t-,o(x) — L,R t T;i{x) when x w 0.7067; Uji t T-.a( x ) = U Rt T;i{x) when x w 0.6573; 

Lt,s-,o( x ) = ir,s ; i(a;) when a; w 0.7969; J/t.SjoOe) = Ut,s-,i(x) when a; w 0.7784; 

£r,S ; o(z) = £r,S;i(z) when x w 0.7916; U R , s - fi (x) = U RiS ;i{x) when x w 0.7737. 

Remark 2.5. We note that piecewise quadratic bounds even tighter than the Lt ± ,t 2 an d Ut x ,t 2 could be 
obtained from Theorem 2.1. The bounds on the ARE given in Corollary 2.3 are derived by appropriately 
rewriting the inequalities g a (0+) < q a (x) < q a (l—) for x € (0, 1) and a e {0, 1}. Of course, one may use 
any finite partition = x < x\ < ■ ■ ■ < x n -\ < x n = 1 of the interval (0, 1) to obtain the corresponding 
piecewise quadratic bounds based on the inequalities q a ( x i-i+) < q a {%) < q a ( x i—) for % € (a;j_i,x,), for 
each i = 1, . . . , n. We state this as another corollary, whose proof will be omitted due to its similarity to 
that of Corollary 2.3. 

Corollary 2.6. Let = xq < x\ < . . . < x n -\ < x n = 1. Then for ARE = AREt 1j t 2 an d q a = qT u T 2 -.a with 
(Ti,T 2 ) e {{R,T),(T, S), (R, S)} and a e {0,1}, one has 

L a ^L^ ARE ^U^U a 

on (0, 1), where 

L a (x) := L TuT2 . a (x) := ARE (a) + ARE' (a) (a; - a) + o Q (x 4 _i+)(x - a) 2 

and 

U a {x) := U Ti .t 2 -A x ) ■= ARE(a) + ARE' (a) (a; - a) + «„(a?i-)(a; - a) 2 

/or x G (Xi-i,Xi), and L := L V Li and {/ := t/ A U\. 

Corollary 2.6 is illustrated by Figure 3; the bounds L and U are based on the partition = xq < x\ < 
x 2 = 1, where x\ — ( x i)t u t 2 is chosen as the mean of the solutions to L = L\ and Uq = U\ (from 
Corollary 2.3), whose approximate values are given in Remark 2.4. That is, (xi)_r,t ~ \ (0.7067 + 0.6573) ~ 
0.6820, {x x ) T ,s ~ ^(0.7969 + 0.7784) w 0.7876 and (x^^s « 1(0.7916 + 0.7737) w 0.7826. 




Fig 3. Illustration of Corollary 2.6, using partition = xo < x\ < X2 = 1, where (xi)r,t ~ 0.682, (xi)t,S ~ 0.788, 
(^Oh.S ~ 0.783; piots are o//t 1: t 2 / or / = ARE (solid), f = L (dotted) and f = U (dashed) 



Note also that Corollary 2.3 immediately implies even better quartic bounds on ARE^s: 
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Corollary 2.7. Let 



Lr,T;o, ■ Lt,S;. 



a 



and U R . s -a ■= U RiT -a ■ U T ,s v 



a 



for a e {0, 1}, and also let 



Lr,s '■= LR t s;0 V L R> s-i and U R> s '■— U R s-,o A U R , s -i- 



Then 



L R ,S;a < L R , S -a < ARE^s < U R ,S;a < U R S-, 



for a € {0, 1} and 



Lr,s < L Rj s < AKE RiS < U Ri s < U Rj s 



on (-1,0) U (0,1). 
3. Proofs 

We first provide Berry-Esseen bounds for the distributions of the test statistics R, S, and T and explicit 
expressions for the asymptotic mean and variance for each of these statistics. Once these facts are established, 
Theorem 2.1 will be proven with the aid of l'Hospital-type rules for determining the monotonicity pattern 
of a ratio. 

3.1. Berry-Esseen bounds and expressions for the asymptotic means and variances of R, S, 
and T 

Each of R, S, and T shall be shown to satisfy (1.1). For each of these statistics, it will be clear that V 
in (1.1) may be taken to be any open interval containing po whose closure does not contain the points 
— 1 or 1. Further note that each of these three statistics is invariant to linear transformations of the form 
Xi i ► aXi + b and Y{ cYi + d with a > and c > 0. So, let us assume without loss of generality (w.l.o.g.) 
that Hx = [i-Y = and ax = oy = 1. For convenience we allow the values p = ±1; then the bivariate normal 
distribution is degenerate: Yi = ±JQ a.s. 

Based on the results of [20, (4.9)], one has the following uniform Berry-Esseen bound on the distribution 



of R: 




here and in what follows, A stands for different positive absolute constants, and 



(3.1) 



a%(p) := E^XM l(Xf + Y?)f = (1 - p 2 ) 2 ; 



the last equality can be checked using the representation Yi = pX\ + \J\ — p 2 Zi, where X\,Z\ ~ iV(0, 1). 
Thus, (1.1) holds for R, with 



(3.2) 



By (2.3), T is a [/-statistic with kernel hr of degree m = 2. Further, S (defined in (2.2)) is a ^-statistic 
with a kernel of degree m — 3; Hoeffding [6, Section 5c] describes how any F-statistic can be expressed as 
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a [/-statistic of the same degree, so that S is a [/-statistic with a symmetric kernel hs, n of degree m = 3. 
Namely, 

S- 

I -) 

i<j< 



V3/ l^i<j<k(^n 



where 
(3.3) 

(3.4) 
and 



n- 2 



VnW, V$, V fc ) := — — y fe ) + — — (MK, V-) + hriVi, V k ) + M^fc))> 

71 | _1_ 77/ | J. 

h s := 2(Kijk + K ikj + Kjik + Kj k i + Kkij + K kji ) - 3 



K ijk := \{Xj < X,} \{Y k < Yi}. 
It follows by Chen and Shao's result [3, (3.4) in Theorem 3.1] that for m € {2, 3} and n> m 



(3.5) 



sup 



n /U-EU 

P < Z 



)-$(*) 



AC 3 



where [/ = (^E i<*i<— <»m^n M^u> • • • ' ^»m) i s an y [/-statistic with a symmetric kernel h such that 
< C for some constant C > 0, 

a\ := m 2 Var#(Vi) > 0, 

and 

ff (Vi) := E[/i(Vi,..., V TO )|Vi]. 
Now consider T as expressed in (2.3), and recall that |/it| = 1. One also has 



(3.6) 



Ht(p) ■= E P T = E p h T {V u V 2 ) = 2 E p ( J 12 + J 21 ) - 1 = 4 E p J 12 - 1 = f sin" 1 p. 



In order to see this, note that E p J 12 = P p {X 1 - X 2 > 0, Y± - Y 2 > 0) = P(Zi > 0, pZi + ^1 - p 2 Z 2 > 0), 
where Z\ and Z 2 arc independent standard normal r.v.'s. By the circular symmetry of the distribution of 
(Zi, Z-i) on the plane, we see E p Ji 2 is simply the proportion of the length of the arc of the unit circle between 
the points (0, 1) and — p 2 , —p); that is, 

EpJi 2 = £(f -sin" 1 (-p)) = i + ^-sin-V, 

whence (3.6). 

One can use a similar geometric reasoning to obtain an expression for the asymptotic variance of T. Let 

0r(Vi) :=E[M^i.^2)|^i] =2E[Ji 2 + J 2 i|Vi] -1, 

so that 

o\(p) := 4 Var p fl T (^i) = 16 (E p J i2 J13 + 2 E p J 12 J 31 + E p J 21 J 31 - 4[E p J 12 ] 2 ). 
Consider first E p Ji 2 J13 = P([/i > 0, 17 2 > 0, U3 > 0, U4 > 0), where the Uj's are standard normal r.v.'s with 



£ := Cov 







"1 


1 
2 


P 


P" 

2 


C^ 2 




1 

2 


1 


P 
2 


P 


^3 




P 


P 
2 


1 


1 
2 


[/ 4 




P 
.2 


P 


1 

2 


1 



Raymond Molzon and Iosif Pinelis/Monotonicity properties of the ARE 



9 



That is, E p J12J13 is the probability that the random point S 1 / 2 [Zi, Z 2 , Z 3 , Z^ lies in the first orthant 
of 4-dimensional space, where the Zj's are independent standard normal r.v.'s; further, this is simply the 
ratio of the volume V(p) of the spherical tetrahedron A\A 2 A 3 A^ to the volume 2-7T 2 of the unit sphere 
S3 := {x € K 4 : ||a;|| = 1}, where the vertices Ai,A 2 ,A 3 ,A 4: of the tetrahedron are the columns of £~ 1//2 
normalized to be unit vectors. One can use the classical result of Schlafli [22] to obtain the volume of 
this spherical tetrahedron. But, in fact, this work has been indirectly done by David and Mallows in their 
derivation of the variance of S; the probabilities E p J\ 2 J\ 3 and E p J12 J31 correspond to correlation matrices 
(r) and (w), respectively, in Appendix 2 of [5]. Using the formulas there, and noting E p J21J31 = E p J\ 2 J\ 3 
by the symmetry of the normal distribution, one sees 

(3-7) 4(/>) = ^-^(sin- 1 f) 2 , 

which is bounded away from over any closed subinterval of (—1,1), so, by (3.5), one has (1.1) for any 

o = poe(-i,i). 

We remark that Kendall's monograph [9, Chapter 10] contains derivations of (3.6) and (3.7). Further, 
Plackett [21] has obtained a more general method for calculating P(Ui > a\,U 2 > 0-2, U 3 > 03, [4 > a 4 ) 
which reduces to the Schlafli method when the cti are all 0. 

Directing attention to S, first note that /ig ; „ is bounded (in fact, one can check that \hs n {Vi, V 2 , V3) | £ 
{1, ^fpj } a.s.). Using geometric reasoning similar to that used to compute E p J i2 (only now using the fact 
that X\ — X 2 and Y\ — Y3 have a correlation of |) , one finds 

EpATi23 = l + ^sin- 1 !, 



so that 

MAP) ■■= E P S = E p ^(^1,^2,^3) = f sin" 1 * -l 3 -^- 

accordingly, let 



2 1 rs+1 



— 6 £ 



(3.8) Ms(p) := lim /x s ,„(p) = E p h s (V u V 2 , V 3 ) = f sin 

n— ►oo 

and note that \fn(psn — Ms) — > uniformly for p e (—1, 1). 
Let next 



9S.n (Vi) 

ffs(Vi) 

<4(p) 



= E[/is,n(Ui, U 2 , U 3 )|Vi] ; 

= E[MVl,V2,V 3 )|Vl] =4E[K 123 + ^213 + ^23l|^l] -3; 

-9Var p5s ,n(Ui); 
= 9Var p < 7s (F 1 ) 

= 144(E p X 123 Xi45 + 2 E p K 213 ^4i5 + 4 E p ^ 123 ^ 4 i 5 + 2 E p X 213 ^45i - 9[E p K 123 } 2 ) 



where (Ts,n(p) ■= \J o"|,„(p) and fT s(p) := noting that E p K 231 K ibl = E p K 213 K 415 and E p K 123 K 451 

= E p K132-K415 = E p K 123 K 4 i 5 since the distributions of ((X 1: Yi), (X n , Y n )) and ((Yi, X x ), (Y n , X n )) 
are identical and permutation-invariant. It is clear that expressions for <r| n and cr| may be derived in terms 
of the volumes of spherical tetrahedra via Schlafli 's formula. For the sake of brevity, we omit these details and 
refer the reader to David and Mallows' derivation of Var S; note the probabilities E p K \ 23 Ki 45 , E p Ki 23 K 4 i 5 , 
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Ep K213K415, E p K213K451 correspond to the correlation matrices (c), (d), (e), and (f), respectively, found in 
Appendix 2 of [5]. Then one has 



(3.9) 



4(p) 



324 



sin 



2 72 , 
+ -2 (h(p) 



2I 2 (p)+2I 3 (p)+4I 4 (p)), 



where 



/,,:,-,: / Sm ^!ZT 2 ) dv. 
Jo 



-1 u 



f 2(3 



3 sm 



m(4-m 2 ) 



2y2V8-6w 2 +tt 4 



sin" 

== / — 
Jo 



-1 m(4-m 2 ) 
2V12-7m 2 +m 4 



an explicit expression of a Sn is not of direct concern to us and so is omitted (though could also be obtained 
from [5]). Note the integrals are expressed differently than the corresponding ones found in [5], 

though a simple change of variables shows their equivalence (the integrals equivalent to Ii\I 2 \I 3 \ I4 are found 
in Appendix 2 of [5] in the expressions corresponding to the correlation matrices labeled there by (f); (c); 
(f); (d) and (e), respectively.) 

Now, o~s,n — ► &s uniformly over all p G [—1, 1] (since, by (3.3), hs n ~ = 0(l/n)). It will be pointed 

out in the last paragraph of part (TSO) of the proof of Theorem 2.1 that <r| > for p e (—1, 1). It is also 
clear from (3.9) that tr| is a continuous function of p, so that the minimum of as over any closed subinterval 
of ( — 1, 1) is strictly positive. Thus, inf pg y <Js,n(p) > for all large enough n, where V is as introduced in 
the beginning of Section 3.1. Referring now to (3.5) (and replacing there U with S, EU with ps,n and u\ 
with <Js,n), one finds that 



sup 

zeR 



g-Ms(p) 
°s(p)/Vn 



< z 



) - *(z) 



a s,n(p)Vn 



whcr 



(z + 



Ms(p)-MS,n(p) ' 



; in turn, the last two terms in the above inequality vanish uniformly 



<rs(p)/\/n 

over z € R and p e V as n tends to 00 (using well-known properties of the function $ and the previously 
noted facts that \/n(ps — Ps.n) — ► and as, n /o~s — * 1 uniformly on V), so that S 1 satisfies (1.1). 
The next result will be used in the proofs of the statements (TSO) - (RSI) in Theorem 2.1: 



Lemma 3.1. One has a s (l- 



0. 



Proof. W.l.o.g., Yi = pXi + a/1 — p 2 Zi for all i, where the Zj's are i.i.d. AT(0, 1) r.v.'s independent of the 
XiS. Further note that a Sn (p) differs only by a positive constant factor from Var p proj £ S, where C is the 
space of all linear statistics. Also, for p = 1, one has 5 = 1 a.s. and hence Var p proj £ S ^ Var p 5 = 0, so 
cr| n (l) = for all n. Now, letting n — ► 00, one has cr|(l) = 0, since hs. n — hs = 0(l/n). 

Next, \a%(p) = Vir p g s {Vi) = E p h s (W u W 2 , W 3 )h s (W u W i} W 5 ) ' - E 2 p h s (W u W 2 , W 3 ), with W t := 
Wi(p) :— (Xi,pXi + \]l — p 2 Zi). Next, hs(W\, W 2 , W3) and W4, W5) are continuous in p on the 

complement of the union of all events of the form {Xi = Xj} for i ^ j. The latter union has zero probability. 
So, by dominated convergence, o~ s (p) —* c s (l) = as p | 1. □ 

While the result of this last lemma should not be surprising, it should be noted that trying to assert 
ct|(1— ) = using only the expression (3.9) is a more difficult task. 
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3.2. Proofs of monotonicity 

As in [17, 18, 19], let — oo ^ a < b ^ oo, and suppose that / and g are differentiable functions on (a, b). Let 
r := ^ and p := j^; from hereon, the symbol p should not be considered the correlation of a bivariate normal 
population, which latter will be denoted by x. Assume that either g < or g > on (a, b), and also that 
g' < or g' > on (a, b). For an arbitrary function ft defined on (a,b), adopt the notation "ft /*" to mean 
ft is (strictly) increasing on (a, b) and similarly let "ft \" mean ft is decreasing on (a,b); the juxtaposition 
of these arrows shall have the obvious meaning, e.g. "ft / \" means that there exists some c e (a, 6) such 
that ft y on (a, c) and ft, \ on (c, b). Further, let the notation "ft is H — " mean that there exists c <5 (a,b) 
such that ft, > on (a, c) and ft < on (c, 6); similar meaning will be given to other such strings composed 
of alternating "+" and "— " symbols. 

Theorem A: Special-case (l'Hospital-type monotonicity) rules. 

Suppose that either f(a+) = g(a+) = or f(b—) = g(b—) = 0. 

(i) If p f on (a, 6), then r' > on (a, 6) and hence r f on (a, 6); 
(ii) If p \ on (a, 6), £/ien r' < on (a, 6) and hence r \ on (a, 6). 

Theorem B: General (l'Hospital-type monotonicity) rules. 

(%) If p /" and gg' > on (a, b), iften r \, r / or r \ /" on (a, 6); 

fiij // /J /" and gg' < on (a, 6), <ften r \, r or r /* \ on (a, b); 

(Hi) If p \ and gg' > on (a, 6), then r \, r /" or r / \ on (a, 6); 

(w,) If p \ and gg' < on (a, &), <ften r \, r orr \ / on (a, 6). 

Theorem C: Refined general (l'Hospital-type monotonicity) rules. 

Let p := 5 2 = sign( 5 ')(p.g - /). 

fij I/p and gg 1 > on (a, 6), £ften p /; 
(ii) If p y and gg' < on (a,b), then p \; 
(Hi) If p \ and gg' > on (a, 6), £ften p \; 
fraj If p \ and gg' < on (a, 6), then p / . 

In addition, sign(p) = sign(r'), so that the monotonicity pattern of r may be determined by the monotonicity 
of p and knowledge of the signs of p~(a+) and/or p(b—). 

E.g. suppose it can be established that p / and gg' > on (a, b); if one also knows that r(a+) = — oo 
then the general rules imply r / '. Alternatively, p / and gg' > imply p /; if it can be established that 
p(a+) > 0, then p > on (a, b) and hence r /" on (a, b). We shall make frequent use of these rules throughout 
the proof of Theorem 2.1. The special-case rules are proved in [18, Proposition 1.1], and a proof of the general 
rules is found in [17, Proposition 1.9]. A proof of the refined general rules, along with several other variants 
of these monotonicity rules, is found in [19, Lemma 2.1]. Note that Anderson ct al. [1, Lemma 2.2] proved a 
variant of the special-case rules, wherein the function /* (or \) whenever p y (or \). 

That (1.2) may be used to express any of the three pairwise ARE's has been justified by the work of the 
previous section. The proofs of the six statements (RTO) - (RSI) in Theorem 2.1 will follow the same general 
method. Fix an arbitrary a € [0, 1], and let 

(3.10) b := b(a) := ARE (a) and c := c(a) := ARE' (a). 

Then 

ARE(x) — 6 — c(x — a) f(x) — bg(x) — c(x — a)g{x) 
^ a (x — a) 2 (x — a) 2 g{x) 
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when / and g are functions chosen so that ARE = ^ . Accordingly, let 



f (*^0 

fo(x) := f(x) - bg{x) - c(x - a)g(x), g (x) := {x - a) 2 g{x), r {x) := = q a {x), 

(3.H) / , i ^ 

fi-=aifi-i, 9i-=<H9i-i, n — —, := -7-^ = rj, and pi = sign(flr i+ i)(r i+ i5i - /i) 



where the are positive on (0, 1). There is some freedom in choosing the functions a,, though the goal is 
to ensure that, for some natural number n > 1, the ratio r n is an algebraic function. In our case it will turn 
out that r n is actually an algebraic function independent of the value of a. As r n is algebraic, the problem of 
determinining its monotonicity pattern on an interval is completely algorithmic (cf. [23, 4]); here, we use the 
Mathematica Reduce command to deduce the monotonicity of r n — p n -\- The specific choices of /, g and 
the di are given in Lemmas 3.2 - 3.4 below. One may refer to this first phase of the proof as the "reduction" 
phase. 

Once the monotonicity of r n = p n -i is established, the second and final stage of the proof is to "work 
backwards" by using the various l'Hospital-type rules stated above to deduce the monotonicity patterns of 
r n -i = Pn-2, r n -2 = Pn-3, ■ ■ ■ , r i = Po, r o = 9a- Throughout the proof, all functions shall be assumed to be 
defined on (0, 1) unless otherwise stated. 

As most of the functions being treated are rather unwieldy, all calculations are performed with the 
Mathematica (v. 5.2 or later) software; detailed output from the notebooks has been reproduced as the 
appendices, and the actual notebooks will be made available upon request. Each of the appendices RT, TS, 
and RS follows the same general format: the first section (labeled RTr, TSr, or RSr - where "r" stands 
for "reduction (phase)") is dedicated to proving the corresponding one of the Lemmas 3.2-3.4 below (i.e., 
the "reduction" stage of the proofs), the second section (RTO, TSO, or RSO) provides numerical support for 
proving the monotonicity of qo, and the third section (RT1, TS1, or RSI) provides support for proving the 
monotonicity of q\. 

We prove q a is increasing only for a € {0, 1}; the following three lemmas could perhaps be used as starting 
points for the "working backwards" phase for other choices of a G (0, 1) to get even more quadratic bounds 
on the ARE's (cf. Corollary 2.3). It is of course desirable to demonstrate that q a / for arbitrary a € [0, 1] 
(should this be true), though a proof of such a statement has yet to be found; for any given a £ (0,1), 
this second phase of the proof is restricted only by computational capacities, since, as mentioned above, the 
expression for r n is eventually algebraic. We remark also that this method could conceivably be adapted (by 
using an appropriate variant of the definition of q a ) to finding quadratic bounds on AREt.r, = 1/ ARE^r, 
ARE,g T = 1/ ARE Ti g, and AREg^ — 1/ ARE^.g, or possibly finding approximating polynomials of degree 
greater than 2. 

Lemma 3.2. Let a e [0, 1] be arbitrary, and let 

f(x) :=7T 2 -36(sin" 1 f) 2 and g(x) := 9(l - x 2 ) 

for x E (0, 1). Then on the interval (0, 1), one has ARE^t = g, ^4 / , fi < and g 4 < 0, where 

r a/4^ (2 -x 2 ) 2 (50 - 29x 2 + 9z 4 ) 2 

a,{x):=^A-x\ a 2 (x) := a 3 (x) := 5Q _ - , a A {x) := — , 



and fi, gt, Ti are as defined in (3.11). 
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Proof. From (3.2) and (3.6), fJ,' R (x) = 1 and fj,' T (x) = f (1 - x 2 )^ 1 / 2 ; then ARE^ T = f - upon recalling 
(1.2), (3.1), and (3.7). Visual inspection shows that a, > on (0, 1) for i — 1, . . . , 4. It is easily verified that 
fi < 0, #4 < and r' 4 > on (0, 1); see Appendix RTr for explicit expressions of these and the intermediate 
functions. □ 

Lemma 3.3. Let a e [0, 1] be arbitrary, and let 

, , , ^ a, ^ , ^ 4(l-a; 2 )(7r 2 -36(sin- 1 f) 2 ) 
(3.12) f{x):=a%{x) and g(x) := ^ ,K _ J ^i, 

w/iere erf is giijen in (3.9). T/ien on i/ie interval (0, 1), one /ias AREt.s = {j> r io / \ fio > and g w > 0, 
where 

f x Aj 2 , \ (4-x 2 ) 5 / 2 (2 + x 2 ) 2 

ai(x) := V4 - x\ a 2 {x) := — — — 5 — , o 3 (a;) := 



(4-x 2 )(38-17x 2 -3x 4 ) : 



and a 4 , . . . ,a 10 are functions rational in x and \/4 — x 2 , w/iic/i are positive and continuous on (0, 1), wi£/i 
/i, gi, and Ti as defined in (3.11). 

Proof. The proof is found in Appendix TSr. □ 
Lemma 3.4. Let a e [0, 1] be arbitrary, and let 

.2^ ^ . 36(1 - x 2 ) 2 



.f(x) := o s {x) and g{x) := 



tt 2 (4-x 2 ) ' 



where <t| is groen in (3.9). T/ien on i/ie interval (0, 1), one /ias ARE^s = r§ / , f$ > 0, and 9/5 > 0, 
where 

oi (a;) := \/A-x 2 , a 2 (x) := (4 - x 2 ) 5/2 , a 3 (x) := 



x(41 - 20x 2 + 3x 4 ) ' 

and a,4, as are rational functions, which are positive and continuous on (0, 1), with fi, g ir and as defined 
in (3.11). 

Proof. The proof is given in Appendix RSr. □ 

Before proving Theorem 2.1, recall the implications of (3.11). If on some open subinterval of (0, 1) one 
has fi > (or /j < 0), then on this subinterval /j_i f (or /j_i \), and similarly for the g^s. If gi has fc 
roots in (0, 1), these shall be denoted by Xij, j = 1, . . . , k, with the assumption that x^i < ■ ■ ■ < Xi^; if gi 
has only a single root in (0, 1), it will simply be denoted by Xj. Similarly, the roots of /j whenever they exist 
will be denoted by yi,i,yi,2, ■ ■ ■ (or simply yi if fi has a single root), and if ever r- is shown to have a root 
in (0,1) (there will only be at most one root in what follows), this root will be denoted by Zi. Numerical 
approximations of any of these roots are not of direct concern to us, but rather their positions relative to 
other roots. Such information is easily obtained from evaluation of the respective functions at specific points; 
for instance, if at some step we deduce that f\ and gi are both + — , with /i(0.5) > > <7i(0.5), then it is 
inferred that X\ < 0.5 < y\ (and further, that r\(x\— ) = ^(x\— ) = oo and r\(xi+) = ^(x\+) = -co). 

Proof of Theorem 2.1, (RT0). See Appendix RT0 for more details of the following arguments. Adopt the 
notation of Lemma 3.2, with a = 0, so that, in accordance with (3.10), 6 = ARE# ; t(0) = \ and c = 
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ARE' R T (0) = 0. Noting that f 3 (0+) = g 3 (0+) = 0, one has f 3 < 0, g 3 < (since, by Lemma 3.2, / 4 < and 
34 < 0), and also, by the special-case rules, p 2 — r 3 / (since, by Lemma 3.2, p 3 = r 4 /). 

Next, g 2 \ (as g 3 < 0) and 32(0+) > > 32(1—) imply g 2 > on (0,2:2) and g 2 < on (x 2 ,l); 
similarly, f 2 \ and /2(0+) > > /2(1— ) imply / 2 > on (0, j/ 2 ) an d /2 < on (2/2,1)- Verifying that 
32(0.41) < < /2(0.41), one has x 2 < 0.41 < y 2 , further implying r 2 (x 2 —) = 00 and r 2 (x 2 +) = —00. Noting 
the sign of g 2 g 2 (which is the sign of 3233) on each of (0, x 2 ) and (x 2 , 1), the general rules imply pi = r 2 / 
on each of these two intervals. 

Next, 31 / \ on (0, 1) (as 32 is +— ) and 31 (0+) = > 31 (1— ) imply the existence of a single root x\, with 
x 2 < x\; similarly, f\ /* \ and /i(0+) = > /i(l— ) imply the existence of a single root 31, with 32 < 3i- 
The special-case rules imply r x / on (0,x 2 ) (as /i(0+) = 31 (0+) = 0). Further, 31 (0.71) < < /i(0.71) 
implies x\ < 31, which in turn shows r\(x\— ) = 00 and ri(xi+) = —00; noting the sign of 313^ on each of 
the intervals (x 2 ,xi) and (x\, 1), the general rules imply r\ /* on these two intervals. The continuity of r\ 
at x 2 implies po — r\ / on (0,a;i) and {x\, 1). 

Finally, / (0+) = 3o(0+) = /o(l-) = 3o(l-) = imply both 30 > on (0,1) (since 31 is H — and 
hence 30 /" \ on (0,1)) and r /" on each of the intervals (0,a;i) and [x\,X) (by the special-case rules); 
the continuity of vq at x\ implies gR.r-.o = r o /'on (0,1). Further, the l'Hospital rule for limits implies 
ro(0+) = r 2 (0+) and r (l-) = ri(l-). □ 

Proof of Theorem 2.1, (RT1). See Appendix RT1 for more details of the following arguments. Adopt the 
notation of Lemma 3.2, with a = 1, so that 

fc = ARE fl;T (l-) = f(l-) = ^ = ^ 

and 

c = ARE' fl;T (l-) = ^(1-) = (1-) = f ''^'%)^-^\ 

which follows by repeated application of the l'Hospital rule for limits after noting /(l— ) = 3(1—) = 0. 

Next, 33(0+) > > 33(1-) and / 3 (0+) > > /s(l-) along with 33 \ and / 3 \ (since / 4 < 
and 34 < by Lemma 3.2) shows that 33 and f 3 each have a single root x 3 and y 3 , respectively. Also, 
33(0.6) < < /s(0.6) shows X3 < 33 and hence r3(x3— ) = 00 and r 3 (a;3+) = —00. Noting the sign of g 3 g' 3 
on each of the intervals (0, x 3 ) and (x 3 , 1), the general rules imply p 2 = r 3 / on these two intervals. 

Next, 32 / \ (as 33 is H — ) and 32(0+) = 32(1-) = imply 32 > 0, whereas f 2 / \ and /2(0+) < = 
/ 2 (1— ) imply f 2 has a single root 32- The special-case rules imply r 2 / on (x 3 , 1); as p 2 / and g 2 g 2 > 
on (0,x 3 ) and p 2 (0+) > 0, the refined general rules imply p 2 > and hence r 2 /" on (0,0:3). Noting that r 2 
is continuous at x 3 , one has pi = r 2 / on (0, 1). 

Next, 3! and /i(l— ) = 3i(l— ) = imply both g x < and p = r 1 / on (0,1); similarly, 30 \ 
and /o(l-) = 3o(l-) = imply g > and q R . T ;i = r / on (0,1). Lastly, r (0+) = ^[°^| and also 
r (l— ) = r 3 (l— ), which follows by the l'Hospital rule for limits. □ 

Proof of Theorem 2.1, (TS0). See Appendix TS0 for more details of the following arguments. Adopt the 
notation of Lemma 3.3, with a = 0, so that b = ARE Ti s(0) = 1 and c = ARE^ s (0) = 0. Now, g g /, 
/ 9 /, and /g(0+) = 39(0+) = imply / 9 > 0, 39 > 0, and p$ — r 9 / (using the results of Lemma 3.3 and 
the special-case rules) on (0, 1). Also, 3s(l— ) < 0, /s(0+) > 0, and p 8 (0+) < imply g s < 0, f$ > 0, and 
P7 = r 8 \ (by the refined general rules) on (0, 1). Further, /V(0+) = 37(0+) = imply > 0, 37 < 0, and 
Pa = r 7 \ (again by the special-case rules) on (0, 1). 
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Next, (76 \ and 56(0+) > > .96(1—) imply the existence of a single root x 6 ; f 6 / and /e(0+) > imply 
fa > on (0, 1). The refined general rules imply r 6 / on (0,x 6 ) (as /?6(0+) > 0), and also that pe \ on 
(x 6 , 1). As x 6 < 0.75 (since 56(0.75) < 0), note that pe(xe+) > ps(0.75) > > /5 6 (1 — ) implies r 6 /\ on 
(x 6 , 1). That is, r' 6 has a single root z 6 , and hence we have p 5 = r 6 / on each of (0, Xe) and (x6, Ze) and \ 
on (z 6 ,l). 

Next, 35(0+) > and 35(1-) > (along with .95 / \) imply 55 > on (0, 1); also, / 5 (0+) > implies 
/5 > on (0, 1). As x 6 > 0.5 (since .96(0-5) > 0) and p 5 / on (0,X6) (by the refined general rules), one has 
Ps(0+) < < (05(0.5) < p${xq+); that is, r 5 X,/ 1 on (0,x 6 ), or r' 5 has a single root 2:5 (with 25 < x 6 ). Recall 
that f 5 , /g and 55 are all positive on (0, 1), and also g' 5 < on (x 6 , 1). Then r' 5 — hs^hfh > q anc j nence 

r 5 on (x 6 , 1). (Let us remark at this point that the l'Hospital-type rules could, in principle, be used to 
establish the monotonicity of r§ on each of (x§,z§) and (ze, 1); however, this would necessitate proving that 
p5(z&) > 0, a task which requires more work than simply requesting the Mathematica program to evaluate 
the function at the approximation of the root z 6 .) As r 5 is continuous on (0, 1), we have p 4 = r 5 \ on (0, 25) 
and /* on (25, 1). 

Next, 54(0+) = —00 < < .94(1—) and /4(0+) = —00 < < ) imply the existence of roots x 4 

and 2/4 (as 55 > and f$ > 0). As 54(0.3) < < r' 5 (0.3), we see that X4 > 0.3 > z 5 ; the refined general 
rules imply p 4 / on (0,2:5), and so, ,04(0+) = implies r 4 / on (0,25). Also, 34(0. 4) > > /4(0.4) implies 
X4 < 0.4 < 2/4, so that ^(x 4 — ) = 00 and r 4 (x4+) = —00. The general rules then imply r± / on each of 
(25, X4) and (X4, 1). Further, the continuity of r 4 at Z5 implies p% = r 4 / on both (0,^4) and (X4, 1). 

Next, 53 \y and 33(0+) = < .93(1—) imply the existence of a single root x 3 ; at that, x 3 > x 4 ; 
similarly, /3(0+) = < /s(l— ) implies the existence of y 3 . The special-case rules imply r 3 / on (0,0:4); 
53(0.64) > > /3(0. 64) implies X3 < 0.64 < 53, or ^(2:3— ) = 00 and ^(23+) = —00, so that the general 
rules show that r?, / on (x4,X3) and (£3, 1). As r-^ is continuous at X4, one has P2 — ^3 / on (0,0:3) and 
(ara.l). 

Next, 32 \y, along with 52(0+) > > 52(0.5) and 52(1 — ) > 0, implies the existence of two roots X2.1 
and 2:2,2; similarly, .^(0+) > > /2(0.5) and /2(1— ) > shows / 2 has two roots 1/2, i>2/2, i- Noting that 
52(0.35) < < / 2 (0.35) and also 52(0.86) > > / 2 (0.86), we have 2:2,1 < 0.35 < y 2 ,i < 0.5 < 2:3,2 < 0.86 < 
52,2, whence ^(2:2,1— ) = riix-ifi—) = 00 and r2(x2,i+) = ?"2(x 2 ,2+) = —00; the general rules then imply that 
T2 / on each of (0,2:2,1), (2:2,1,2:3), (2:3,2:2,2) and (2:2,2, 1). The continuity of r 2 at 2:3 implies p\=r 2 / on 
(0,X2, 1), (x 2 , 1,2:2,2) and (x 2 , 2 ,l). 

Next, /i(0+) = 5i(0+) = ) = 5i(l— ) = (together with / 2 and 52 both H h) implies the existence 

of roots Xi and y\. That r\ / on (0,x 2 ,i) and (x2,2,l) is implied by the special-case rules; that r\ / on 
(£2,1, xi) and (xi,X2,2) is implied by the general rules upon noting that 51(0. 62) < < /i(0.62) (and hence 
xi < 51, or ri(xi— ) = 00 and ri(xi+) = —00). The continuity of r\ at X2,i and X2,2 implies po = r i / ° n 
(0, xi) and (xi, 1). 

Lastly, ./o(0+) = 5o(0+) = /o(l — ) = 5o(l — ) = shows 50 > on (0, 1) and also, by the special-case rules, 
r / on (0, Xi) and (xi, 1). The continuity of r at x 4 shows qr.s-.o — r o /"on (0, 1). Further, the l'Hospital 
rule for limits yields ro(0+) = r 2 (0+) and ro(l— ) = r 2 (l— ). 

As promised in the remarks preceding Lemma 3.1, we show that as > on (0, 1) (and hence on (—1,0) 
as as is even). Note /o > (as /o /"\ and /o(0+) = /o(l— ) = 0); by (3.12) and (3.11), and recalling that 
6=1 and c = 0, one has /o = a 2 s — 5, so that a s > 5 on (0, 1). As x 2 g(x) = go(x) > 0, it follows that a 2 s > 0. 
Further note that there is no circular reasoning here; the above proof stands on its own, regardless of any 
probabilistic interpretation we give to the functions / or 5. □ 

Proof of Theorem 2.1, (TS1). See Appendix TS1 for more details of the following arguments. Adopt the 
notation of Lemma 3.3, with a = 1, so that /(l — ) = .9(1—) = /'(! — ) = 5'(1— ) = 0, and repeated application 
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of the l'Hospital rule for limits imply 

(3.13) 6 = ABEt,s(1-) = i (1-) = £(1-) = ftH = 997I k 135 ^ 

and 

r — ART?' n ^ - f'g-fg' n 1-/V/j!(i \ _ f"'g+f"g'-f'g"-fg"' n \ 

C - A «^T,sU J - ^ U J - 2gg> I 1 ") ~ 2( 9 ') 2 +2ff S " ^ J 

1 j _ /W fl +2/"' fl '-2/' g'"-fgW (i \ _ /"'(l-)fl"(l-)-/"(l-)fl"'(l-) 

6g>g"+2gg>" V 1 > 3g"(l-) 2 

Then /g(0+) = 59(0+) = (and /10 > 0, 510 > 0, by Lemma 3.3) imply that fa > 0, 39 > and p 8 = r 9 /" 
(by the special-case rules). Also, /s(0+) > 0, 38(1—) < and ps(0+) < imply fa > 0, g$ < 0, and (by the 
refined general rules) p 7 — r 8 \ on (0, 1). 

Next, 37(0+) > > .97(1—) implies the existence of a single root 2:7; fa(0+) > shows that fa > 0. The 
refined general rules imply p-j /* on (0, x 7 ) and \ on (x 7 , 1). As p 7 (0+) > 0, we see r 7 f on (0, x 7 ); further, 
x 7 < 0.2 (implied by 37(0. 2) < 0) yields /5 7 (x 7 +) > p 7 (0.2) > > /5 7 (1-), so that r 7 /" \ on (x 7 , 1). That 
is, pe — r 7 / on both of (0,x 7 ) and (x 7 ,z 7 ), and ps = r? \ on (zj, 1). 

Next, 36(0+) > > 36(1—) implies the existence of xq; /e(0+) > implies fa > on (0, 1). As pe(0+) > 0, 
the refined general rules imply r§ /* on (0,x 7 ). Further, 36(0.5) > > r 7 (0.5) implies z 7 < 0.5 < x$; as 
fa > 0, /g > 0, 36 > 0, and g' 6 < on the interval (x 7 ,x 6 ), we have r' 6 = J(i96 ~ 2 > and hence r e / 

on (xj,xq), so that r 6 / on (0,x 6 ) (since r 6 is continuous at x 7 ). Also, p6 \ on (x 6 ,l) is implied by the 
refined general rules; then 3 6 (0.85) < implies x 6 < 0.85, so that pe(xa+) > p 6 (0.85) > > /5 6 (1 — ) shows 
that r 6 y\ on (x 6 , 1). That is, /O5 = r 6 /"on (0, X6) and (x 6 , z 6 ) and \ on (z 6 , 1). 

Next, 35(0+) > and g 5 (l-) > 0, along with 35 / \, imply 35 > on (0, 1); also, / 5 (0+) < < fa(l-) 
implies fa has a single root 35. The refined general rules imply r 5 /* on (0,xe), as ^5(0+) > 0; also, 
/5(0.5) > implies 35 < 0.5 < x&, so that fa > 0, f 5 > 0, 35 > and g' 5 < on (x 6 ,l), and hence 
r i _ fsgs-hg*. -> on (x 6 , 1)- As r 5 is continuous at x 6 , one has p 4 = r 5 /" on (0, 1). 

Next, —00 = 34(0+) < < 34(1—) shows 34 has a single root x 4 ; /4(0+) = 00 > > /4(0.75) and 
A(l-) > shows / 4 has two roots 34,1 and y 4i2 . Also, 34(0. 75) < < 34(0.8), / 4 (0.75) < 0, and / 4 (0.8) < 
together imply X4 € (0.75,0.8) C (34,1,34.2), so that r 4 (x4— ) = 00 and r 4 (x4+) = —00. The general rules 
then imply p 3 = r 4 / on each of (0, X4) and (x 4 , 1). 

Next, 33(0+) > = 33(1—) and 33 \/* shows 33 has a single root X3; fa(0+) > = /s(l— ) and 
fa / \ / shows fa has a single root 33. Then r 3 / on (x 4 , 1) by the special-case rules; 33(0.5) < < /s(0.5) 
yields x 3 < 33 (and hence r 3 (x3~ ) = 00 and r3(x3+) = —00), so that the general rules imply r 3 / on both 
of (0,2:3) and (2:3, £4). As r 3 is continuous at x 4 , p 2 = r 3 / on (0,2:3) and (x 3 , 1). 

Next, 32(0+) < = 32(1—) and /2(0+) < = fa(\— ) together yield the existence of roots x 2 and 32, 
along with r 2 / on (2:3,1) (via the special-case rules). Also, 32(0.1) > > /2(0.1) implies x 2 < 32 (and 
hence ^(2:2— ) = 00 and ^(2:2+) = —00), so that the general rules then imply r 2 f on (0,2:2) and (2:2,2:3). 
Further, r 2 is continuous at X3 and hence p\—r 2 / on (0,2:2) and (2:2, 1). 

Next, 31(0+) < = 3i(l-) and /i(0+) < = /i(l-) show that 31 < and fa < on (0, 1), and also 
r\ /* on (x 2 , 1) by the special-case rules; pi(0+) > implies via the refined general rules that r\ f on 
(0,2:2). The continuity of n at x 2 then shows p = n / on (0, 1). 

Lastly, /o(l— ) = 3o(l — ) = shows that 30 > and further, via the special-case rules, that Qt,S;1 — r o / 
on (0, 1). Note r (0+) = ^±1 anc j 5 by the l'Hospital rule for limits, r (l-) = r 4 (l-). □ 

Proof of Theorem 2.1, (RS0). See Appendix RS0 for more details of the following arguments. Set a = in 
the notation of Lemma 3.4, so that, in accordance with (3.10), b = ARE^, s (0) = ^ and c = ARE' fi s (0) = 0. 
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Then /4(0+) = 34(0+) = 0, /g > 0, and 35 > (from Lemma 3.4) together imply that f± > 0, 34 > 0, and 
also p 3 — r4 / (via the special-case rules). 

Next, g 3 /* and 33(0+) < < .93(1—) implies the existence of the root X3; that f 3 has a single root 2/3 
follows by f 3 / and / 3 (0+) < < / 3 (1-). From x 3 < y 3 (implied by g 3 (0.64) > > / 3 (0.64)) follows 
^3(^3—) = 00 and ^3(^3+) = — °°; the general rules then imply p 2 = r 3 / on both (0,x 3 ) and (x 3 , 1). 

That gi has two distinct roots X2,i and X2,2 follows from 32(0+) > > 32(0.5) and 32(1—) > (along 
with g 2 \y); similarly f 2 has two roots 2/2,1 and 2/2,2, which follows from /2(0+) > > / 2 (0.5) and 
.Ml-) > 0. Then #2(0.33) < < / 2 (0.33) shows that x 2 ,i < 2/2,1, and g 2 (0.S6) > > / 2 (0.86) (together 
with > (/2(0.5) and > ^2(0. 5)) show that 2/2.1 < 0.5 < X2.2 < J/2, 2- The general rules then imply (since 
^2(22,1-) = ^2(^2,2-) = 00 and r 2 (x 2 ,i+) = ?"2(x2,2+) = -00) that p\ = r 2 / on the four intervals (0,x 2; i), 
(x2,i,x 3 ), (x 3 ,x 2i2 ) and (x 2j2 , 1); the continuity of r 2 at x 3 implies p\=r 2 / on (x 2l i,X2,2). 

As /i(0+) = <7i(0+) = ) = .gi (1— ) = 0, one finds the existence of roots x\ and 2/1 (since g 2 

and f 2 are both H h), as well as r\ / on (0,X2,i) and (£2,2,1) via the special-case rules. Further, 

2/i(0.6) < < /i(0.6) shows x\ < 2/1 (and hence r\{x\—) — 00 and ri(xi+) = —00), so that the general rules 
imply ri / on (x 2 ,i,xi) and (xi,X2, 2 )- The continuity of r\ at x 2j i and x 2 .2 then implies p — r\ / on 
(0, x\) and (xi, 1). 

Lastly /o(0+) = .90 (0+) = /o(l — ) = .9o(l — ) = and go S \ imply that 2/0 > and also (by the special- 
case rules) that ro /* on both (0,Xi) and (x\, 1). The continuity of ro at x\ then implies qR,s-o = ^0 / on 
(0, 1). The l'Hospital rule for limits implies ro(0+) = r2(0+) and ro(l— ) = r 2 (l— ). □ 

Proof of Theorem 2.1, (RSI). See Appendix RSI for more details of the following arguments. Adopt the 
notation of Lemma 3.4, with a = 1, so that /(l— ) = .9(1—) = /'(l— ) = </(l — ) = and repeated application 
of the l'Hospital rule for limits together yield (similar to (3.13) and (3.14)) 

& = ARE fliS (l-)= ^ = 3 -^=^ and c = ARE^(l-) = H'W-rHAH . 

From 2/4(0+) < < 274(1—) and 2/5 > follows the existence of £4; similarly, /4(0+) < < ) and 

/ 5 > imply the existence of 2/4- Then .94(0.8) > > .A(0.8) shows X4 < 0.8 < 2/4, or hence r 4 (x 4 — ) = 00 
and r4(x4+) = —00, and so the general rules imply p 3 — r^ /" on both (0,0:4) and (#4, 1). 

Next, 273 \y (as 274 is — h) and 273(0+) = 00 > — g 3 (l—) yield the existence of x 3 ; that f 3 has a single 
root 2/3 also follows by f 3 \/* and / 3 (0+) = 00 > = / 3 (1— ). The special-case rules imply r 3 / on (x 4 , 1); 
also x 3 < 2/3 follows from 273(0.5) < < /3(0.5) (whence ^(2:3— ) = 00 and ^(2:3+) = —00), and so, the 
general rules imply r 3 / on both of (0,x 3 ) and (x 3 ,x 4 ). Also, r 3 is continuous at X4 and hence p 2 — r 3 / 
on (0, x 3 ) and (x 3 , 1). 

As 272(0+) < = .92(1 — ) and /2(0+) < = / 2 (1— ) (and 273 and f 3 are both H — ), there exist roots x 2 
and 2/2; the special-case rules imply r 2 / on (x 3 ,l). Further, (72(0.1) > > /2(0.1) shows x 2 < 0.1 < 2/2 
and hence r 2 (x 2 — ) = 00 and r 2 (x2+) = —00. The general rules then imply r 2 / on (0,x 2 ) and (x 2 ,x 3 ); the 
continuity of r2 at x 3 then implies p\ = r 2 /* on (0,^2) and (x 2 , 1)- 

One finds that . 9l < and ft < on (0, 1), as <?i(0+) < = ffi(l-) (with c?i \/) and /i(0+) < = 
/i(l— ) (with fx \y), which further imply by the special-case rules that r\ / on (x 2 , 1). Also, pi(0+) > 
implies via the refined general rules that pi > 0, or r\ / , on (0, x 2 ); as r\ is continuous on (0, 1), one sees 
p = n / on (0,1). 

Lastly, /o(l— ) = 5o(l — ) = imply in the first place that 27o > (as 270 \), and in the second place that 
1r„S:\ — r o / on (0, 1) (via the special-case rules). The l'Hospital rule for limits implies r (l— ) = r 4 (l— ), 
and so (0+)>0 implies r (0+) = ^±l. □ 
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Proof of Corollary 2.2. As the ARE's are even functions here, one has ARE'(O) = 0, and hence ARE(a;) = 
ARE(O) + x 2 q (x) for x G (0, 1). Theorem 2.1 shows q / and q (0+) > 0, which imply q > on (0, 1); 
hence ARE /* on (0, 1) as well. The values ARE(0+) and ARE(1— ) are exactly those values of b given at 
the beginning of the proof of each of the six parts of Theorem 2.1. □ 

Proof of Corollary 2.3. The result immediately follows from Theorem 2.1: 

(x - a) 2 q a (0+) < ARE(x) - ARE(a) - ARE' (a) (a: - a) < (x - a) 2 q a (l-) 
=> L a (x) < ARE (a;) < U a {x) 

for all x € (0,1) and a £ {0,1}. Replacing u x" with "— x" in the above inequality when x e (—1,0) and 
recalling the ARE is even yields the desired results. □ 

Proof of Corollary 2.7. Noting that ARE fl , s = ARE fljT • ARE TiS , one has L RiT ■ L T ,s < ARE fl:S < U RT • 
U T ,s on (-1,1) \ {0}. That L RS > L RS and U RS < U RS is easily verifed by noting L RS - L RS and 
U R s — U R s have no roots on (—1, 1) \ {0} and verifying their appropriate signs. □ 
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Appendix RT 

Raymond Molzon and losif Pinelis 

Michigan Technological University 

RTr: Proof of Lemma 3.2 (Reduction of qRj ;a to an algebraic form) 

Assume by default all functions are defined on (0, 1): 
$Assumptions = < x < 1; 

The asymptotic mean and variance of Pearson's R are: 

= x; (* (3-2) in paper *) 
(ct 2 ) r [x_] = (1 -x 2 ) 2 ; (* (3.1) in paper *) 

The asymptotic mean and variance of Kendall's T are: 

2 

Mt[x_] = — ArcSin[x]; (* (3.6) in paper *) 

7T 

4 16 x 2 

(<7 2 ).[x_] = ArcSinf-1 ; (* (3.7) in paper *) 

T 9 n 2 2 J 

Express AREr j(x) as the ratio ARE/; j = ~- 

X 2 

f [x_] = 7r 2 - 36 ArcSinf— 1 ; 

2 

g[x_] =9 (l-x 2 ); 

(° 2 ) T M ^ R .[x] 2 

ARE[x_] = ; 

(a 2 ) R [x] u T ' [x] 2 

f [x] 

ARE [x] == // Simplify 

g[x] 

True 

This confirms that ARE R T = — , 



Begin the reduction phase, letting b = ARE R 7 e = AREJ, r (d) and qnj;a( x ) — - >~o(x), for arbitrary a e [0, 1]. 

Follow the scheme outlined in the paper (cf. (3.1 1) there), wherein / = a, f l '_ l , gj = a, g' j _ l , r, = — = = p,_i, and ci, > on (0, 1). 

St 8j-t 

At each step, note the continuity of fi and g, in x e (0, 1). 

f0[x_] = f [x] -bg[x] - c (x - a) g[x] // Simplify 

7T 2 + 9b(-l + x 2 )-9c(a-x) (-1 + x 2 )- 3 6 ArcSin [ — ] 

gO[x_] = (x-a) 2 g[x] //Simplify 

-9 (a - x) 2 (-1 + x 2 ) 



1 
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fO[x] 

rO[x_] = ; 

gO[x] 

Note that b and c are left as unspecified constants in the reduction phase; they will be given specific values (upon a choice of a) in the "working backwards" phase 
of sections RTO and RT1 . 



al [x_] = ^4 - x 2 ; 

This choice of a\ yields a constant coefficient on sin~'( in f\. 



fl[x_] = Map[Simplify, Collect [al [x] f ' [x] , {ArcSin[_], Sqrt[_]}], 2] 

-9 \Ja - x 2 (c-2bx + 2acx-3cx 2 ) -72 ArcSin [ — j 
gl [x_] = al [x] gO ' [x] // Simplify 



-18 \Ia - x 2 (a - x + a 2 x - 3 a x 2 + 2 x 3 ) 
fl[x] 



rl [x_] 



gl[x] 



A/4 -x 2 
a2 [x_] - -J 



2 -x 2 

This choice of a-i yields a constant coefficient on b in fa. 

f2[x_] = Map[Simplify, Collect [a2 [x] f 1 ' [x] , {b, c, a}], 2] 

72 9 c (8 a - 25 x - 4 a x 2 + 9 x 3 ) 

36 b + + 

-2 + x 2 -2 + x 2 

g2[x_] = Map[Simplify, Collect [a2 [x] gl ' [x] , a], 2] 

18 a x (-25 + 9 x 2 ) 36 (2 - 13 x 2 + 4 x 4 ) 



-36 a 2 



r2 [x_] 



-2 + x 2 -2 + x 2 

f2[x] 



g2[x] 



Note that ri is a rational function, and indeed one could begin using common analytical techniques to deduce the roots and monotonicity patterns of /?, g?> and i~2- 
We continue with the reduction phase, so as to obtain a ratio independent of the value of a (and hence b and c); this is done primarily for the sake of consistency 
with the reduction phase in Appendices TS and RS (wherein more complicated functions are considered). 

(2-x 2 ) 2 

a3 [x_] 



50 - 29 x 2 + 9 x 4 

This choice of yields a constant coefficient on c in fa and a in g3 : 

f3[x_] = Map[Simplify, Collect [a3 [x] f2 ' [x] , c] , 2] 

144 x 

9 c 

50 - 29 x 2 + 9 x 4 

g3[x_] = Map[Simplify, Collect [a3 [x] g2 ' [x] , a], 2] 

288x (6-4 x 2 +x 4 ) 



18 a 



r3 [x_] 

g3[x] 



50 - 29 x 2 + 9 x 4 
f3[x] 



2 
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(50 - 29 x 2 + 9x 4 ) 2 

a4[x_] = ; 

2 -x 2 

f 4 [x_] = a4 [x] f 3 ' [x] // Simplify 

-144 (25 + 27 x 2 ) 

g4 [x_] = a 4 [x] g3 ' [x] // Simplify 

8 64 (-50 + 4 6 x 2 - 11 x 4 + 3 x 6 ) 

f4[x] 

r4[x_] = ; 

g4[x] 

Now f 4 , g 4 and r 4 are all algebraic functions, independent of the value of a. 

Reduce [f 4 [x] < SS < x < 1] 

< x < 1 

Reduce [g4[x] <0 && < x < 1] 

< x < 1 

Reduce [r4 ' [x] > && < x < 1] 

< x < 1 

Thus, f\ < 0, g4 < and r\ > on (0, 1), so that = pi = — /* on (0, 1). Note this is true for arbitrary a e [0, 1]. This proves Lemma 3.2 in the paper. 



3 
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RTO: Monotonicity of qRj ; o on (0, 1) 

Below are numerical calculations supporting the arguments of the proof of statement (RTO) in Theorem 2.1. 
{a, b, c} = {0, ARE [0] , ARE ' [0] } 

{o, — , 0} 

L 9 J 

{g3[0], f3[0] }// Simplify 

{0, 0} 



Sign /@ {g2[0] , g2[41/100], g2[l]} 

{1, -1, -1} 

Sign I® {f2[0] , f 2 [41 / 100] , f2[l]} 

{1, 1, -1} 

Sign l@ {gl[0] , gl[71/100], gl[l]} 

{0, -1, -1} 

Sign I® {fl[0] , fl[71/100], fl[l]} 

{0, 1, -1} 

{g0[0], g0[l]} 

{0, 0} 

{f0[0], f0[l]} 

{0, 0} 

{r2[0] //Simplify, r2[0] / / N} 

r 1 i 
|- (-9 +7T 2 ) , 0.0966227j 

{rl[l] //Simplify, rl[l] / / N} 
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RT1 : Monotonicity of «7/?,r ; i on (0, 1) 



Below are numerical calculations supporting the arguments of the proof of statement (RT1) in Theorem 2. 1 . 

{f [1], g[l]} 

{0, 0} 

- f'[l] f "[1] g'[l] -f '[1] g"[l], 

{a, b, c} = 1, , //Simplify 

g'[l] 2g'[l] 2 1 

fl, -— (-9 + V3~^)} 

1 3-v/T 271 ,S 



Sign I® {g3[0] , g3[3/5], g3[l]} 

{1, -1, -1} 

Sign /@ {f3[0] , f3[3/5], f3[l]} 

{1, 1, -1} 



{g2[0], g2[l]} 

{0, 0} 

Sign /@ {f2 [0] , f2[l] //Simplify} 

{-1, 0} 

Sign[r3[0] g2[0] -f2[0]] (* p 2 (0 + )>0 *) 

1 



{gl[l], fl[l]} 

{0, 0} 



{gO[l], fO[l]} 

{0, 0} 



{r0[0] //Simplify, r0[0] / / N} 

( (l8 - 8 -/Tit + 3 jt 2 ) , 0.1510 

2 7 ' ' 

{r3[l] //Simplify, r3[l] / / N} 

1 



+ 5 V3 7T , 0.224778 
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Appendix TS 
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Michigan Technological University 

TSr: Proof of Lemma 3.3 (Reduction of qr,s;a( x ) to an algebraic form) 

By default, all functions shall be assumed to be on (0, 1): 

$Assumpt±ons = < x < 1 

< x < 1 



in _1 |y3/8 1 Lettingff = sin-'(l/4)and/3 = sin-'|y3/8 1 



A term that shall recur several times is sin (l/4) + 2sin -W3/8 . Letting a = sin (1/4) and /? = sin , note that 



sin (a + 2 /3) = sin (a) cos (2 fi) + cos (a) sin (2 (i) 

= sin (a) (l - 2 sin 2 (13)) +2 cos (a) cos (/3) sin (/3) 
1 15 

16 16 

Thus, a + 2 j3 = — + Ink for some natural number k; that k = is easily numerically verified : 

ArcSin[l / 4] + 2 ArcSin[y3 / 8 ] / / N 

1 . 5708 

7T/2 / / N 

1 . 5708 



sub = |ArcSin[l/4] -» 2 ArcSin ["^3 / 8 ] j ; 



The asymptotic mean and variance of Kendall's T are: 
2 

Mt[x_] = — ArcSin [x]; (* (3.6) in the paper *) 



4 16 r x.^ 

(a 2 ) [x_] = ArcSin — ; (* (3.7) in the paper *) 

T 9 n 2 L 2 J 



1 
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The asymptotic mean and variance of Spearman's S are: 
6 -x, 

H s [x_] =— ArcSin — ; (* (3.8) in the paper *) 
7T L 2 1 

324 r x 2 72 



(ct 2 ) s [x_]:=1 ArcSin[— ] + — (II [x] + 2 12 [x] + 2 13 [x] + 4 14 [x] ) ; (* (3.9) in the paper *) 



ArcSin [ — — 1 ArcSin [ — - — 1 

L 8-4 X 2 J L 6-2 x 2 J 

II ' [x_] = — ^ ; 12 ' [x_] = — ^ ; 

-^4 - x 2 -^4 - x 2 

ArcSin I ' ' I ArcSin I — ' ' I 

2 V2~ *j8-6 x 2 +x 4 2 -W12-7 x 2 tx 4 

13 ' [x_] = — = £ = ; 14 ' [x_] 



^4 - x 2 a/ 



,4 - x 2 

Il[x_?NumberQ] : = NIntegrate [II ' [u] , {u, 0, x}] 
12 [x_?NumberQ] : = NIntegrate [12 ' [u] , {u, 0, x}] 
I3[x_?NumberQ] : = NIntegrate [13 ' [u] , {u, 0, x}] 
14 [x_?NumberQ] : = NIntegrate [14 ' [u] , {u, 0, x}] 
I1[0] =I2[0] =I3[0] =I4[0] =0; 

(a 2 ) g [l] = 0; (* This follows from Lemma 3.1 in the paper. *) 

Express AREj- « as a ratio ARE j s = ~ '■ 

g 

f [x_] := (ct 2 ) s [x] 

4(l-x 2 ) (tt 2 - 36 ArcSin [ |] 2 ) 



g[x_] 



ARE[x_] = 



7T 2 (4-x 2 ) 

(° 2 ) S W ^'[x] 2 
(a 2 ) T [x] Ms ' [x] 2 ' 



f [x] 

ARE [x] == // Simplify 

g[x] 

True 



f 

This confirms that AREj j = - . 



Begin the reduction phase (cf. Appendix RTr and (3.1 1) in the paper), letting b = ARErxfa), e = ARE^ s (a) and qrsix', a ) - ~T7 — foWi f° r arbitrary a e [0, 1]. 



f0[x_] := f [x] -bg[x] - c (x - a) g[x] // Simplify 
g0[x_] = (x-a) 2 g[x] //Simplify; 
f0[x] 

r0 [x_] 



g0[x] 



[x_] = ^4 



al [> 

This choice of a\ ensures all of the sin~'( ■) terms associated with the integrals ly, I 2 , h, and I 4 (found in f ) have constant coefficients. 

fl[x_] = Map[Simplify, Collect [al [x] fO ' [x] //Together, ArcSin[_]], 2] 

gl[x_] = Map[Simplify, Collect [al [x] gO ' [x] //Together, ArcSin[_]], 2] 

fl[x] 

rl[x_] = ; 

gl[x] 
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Some terms that will be convenient to collect in the functions defined from hereon: 



si = -^16 - 12 x' 
s2 = -^16 - 12 x : 



+ x 4 "^12 - 8 x 2 + x 4 
+ x 4 J: 



4 x 2 ; 



s3 = "^16 - 12 x 2 + x 4 "^9 - 4 x 2 ; 
s4 = - Ax 2 -\jl2 - 8 x 2 + x 4 ; 

Below are some terms that will arise in the computation of/?; they will be convenient to rewrite in terms of si, S2, S3, and £4: 



tl = - / (-4 + x 2 ) (-1728 + 3216 x 2 - 2204 x 4 + 676 x 6 - 89 x 8 + 4 x 10 ) 



== (4 - x 2 ) 3/2 si // Simplify 



True 



t2 = A /- ("4 + x 2 ) (192 - 272 x 2 + 124 x 4 - 20 x 6 + x 8 ) 



== (4 -x 2 ) 3/2 s2 // Simplify 



True 



t3 = - / (-4 + x 2 ) (-144 + 172 x 2 - 57 x 4 + 4 x 6 ) 



True 



t4 = - / (-4 + x 2 ) (-108 + 120 x 2 - 41 x 4 + 4 x 6 ) 



(4 - x 2 ) 3/2 s3 // Simplify 



(4 - x 2 ) 3/2 s4 // Simplify 



True 

Then define the appropriate substitutions: 



sub2 = {tl-> (4 -x 2 ) 3/2 si, t2^ (4-x 2 ) 3/2 s2, t3 -> (4-x 2 ) 3/2 S 3, t4 -> (4 - x 2 ) 3/2 34}; 



a2 [x_] 



(4 -x 2 ) ! 



2 + x^ 



This choice of aj yields a constant coefficient on the term b (sin -J in the function j\. 
f2[x_] = Map[Simplify, 

Collect [Map [Simplify, Collect [a2 [x] fl ' [x] // Simplify // Together , {ArcSin[_], si, s2, s3, s4}], 2] /. sub2, 
{ArcSin[_], b, c, si, s2, s3, s4}], 2] 

g2[x_] = Map[Simplify, Collect [a2 [x] gl ' [x] , ArcSin[_]], 2] 
f2[x] 



r2 [x_] 



g2[x] 



(2 + x 2 ) ; 



a3 [x_] 

(4 - x 2 ) h3[x] 

It will start becoming less obvious that the a,'s are positive on (0, 1), and so this will be routinely checked from hereon: 

Reduce [a3 [x] > && < x < 1] 

< x < 1 



This choice of 03 yields a constant coefficient on c (sin -j in/3. 
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f3[x_] = Map[Simplify, Collect [a3 [x] f2 ' [x] //Together, {ArcSin[_], b, c, si, s2, s3, s4}], 2] 

g3[x_] = Map [Simplify , Collect [a3 [x] g2 1 [x] , ArcSin[_]], 2] 

f3[x] 

r3[x_] = ; 

g3[x] 



Some further simplifying substitutions for fa : 



sub2 = Union 



[sub2, | 



-4 + 



\-9 + 4 x 2 



(4-x 2 ) 



2 ,5/2 



(9 - 4x 2 ) f 



h4 [x_] = 892 - 440 x 2 + 61 x 4 - 9 x 6 

(4 -x 2 ) 5/2 h3[x] 2 

a4[x_] = ; 

x (2 + x 2 ) h4 [x] 

Reduce [a4 [x] > SS < x < 1] 

< x < 1 



4 - x^ 



16 - 12 x 2 + x 4 



(4-x 2 ) 



2 v 5/2 



(16 - 12 x 2 + x 4 ) ! 



This choice of a\ yields a constant coefficient on b sin (-J in fa. 

f4[x_] = Map[Simplify, Collect [a4 [x] f3 ' [x] //Together, {ArcSin[_], b, c, si, s2, s3, s4}], 2] /. sub2 

g4[x_] =Map[Factor, Map[Simplify, Collect [a4 [x] g3 ' [x] , ArcSin[_]], 2], 2] 

f4[x] 

r4[x_] = ; 

g4[x] 



h5 [x_] = 328 256 - 60 276 x 2 - 28 380 x 4 + 12 853 x 6 - 678 x 8 + 81 x 10 ; 
x 2 h4[x] 2 



a5 [x_] 



(4 - x 2 ) h3[x] h5[x] 
Reduce [a5[x] > SS < x < 1] 

< x < 1 



This choice of as yields a constant coefficient on the term c sin ( -) in fa . 

f5[x_] = Map[Simplify, Collect [a5 [x] f4 ' [x] //Together, {ArcSin[_], b, c, si, s2, s3, s4}], 2] 
g5[x_] =Map[Factor, Map[Simplify, Collect [a5 [x] g4 ' [x] , ArcSin[_]], 2]] 
f5[x] 



r5 [x_] 



g5[x] 



h6 [x_] = 17 418 976 - 12 356 932 x^ + 3 290 736 x 4 - 575 137 x" + 35 Oil X s - 447 x i0 + 81 x 1 
a6 [x_] 



4 



4-x z h5 [x] ' 



xh4[x] h6[x] 
Reduce [a6[x] > SS < x < 1] 

< x < 1 

This choice of «6 yields a constant coefficient on b in fa; at this point, fa is an algebraic function. 
f6[x_] = Map[Simplify, 

Collect [Map [Simplify, Collect [a6 [x] f5 ' [x] //Together/. sub2, {b, c, Sqrt[_]}], 2] /. sub2 // Together, 
{b, c, si, s2, s3, s4}] , 2] 

g6[x_] =Map[Factor, Map[Simplify, Collect [a6 [x] g5 ' [x] , ArcSin[_]], 2], 2] 
f6[x] 



r6 [x_] 



g6[x] 
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hi [x_] = 

18 745 083 424 - 14 666 397 812 x 2 + 4 272 900 412 x 4 - 473 552 785 x 6 + 47 852 540 x 8 - 89 482 x 10 + 1296 x 12 - 729 x 14 ; 

h6[x] 2 

a7[x_] = ; 

h5[x] h7[x] 

Reduce [a7[x] > SS < x < 1] 

< x < 1 



This choice of a 7 yields a constant coefficient on c in the function fi. 

f7[x_] = Map[Simplify, Collect [a7 [x] f 6 ' [x] //Together, {c, si, s2, s3, s4}], 2] 

g7[x_] = Map[Simplify, Collect [a7 [x] g6 ' [x] //Together, {ArcSin[_], Sqrt[_]}], 2] 

f7[x] 

r7[x_] = ; 

g7[x] 



h8 [x_] = 67 393 220 864 - 66 665 518 536 x 2 + 25 281 966 744 x 4 - 

4 783 210 446 x 6 + 320 370 996 x 8 - 26 281 941 x 10 - 170 777 x 12 + 231 x 14 + 81 x 16 ; 

h7[x] 2 

a8[x_] = ; 

h6[x] h8[x] 

Reduce [a8[x] > && < x < 1] 

< x < 1 



This choice of a 8 yields a constant coefficient on sin 1 ( in the function gg ; note also that all dependence on the value of a has vanished in / 8 . 

f8[x_] = Map[Simplify, Collect [a8 [x] f7 ' [x] //Together, {si, s2, s3, s4}], 2] 

g8[x_] = Map[Simplify, Collect [a8 [x] g7 1 [x] , ArcSin[_]], 2] 

f8[x] 

r8[x_] = ; 

g8[x] 



h9 [x_] = 32 482 389 470 208 - 34 864 017 237 408 x 2 + 16 286 313 144 464 x 4 - 4 430 399 397 672 x 6 + 832 485 830 428 x 8 - 
10 4 5 7 8 2 6 7 9 6 x 10 + 7 8 5 5 4 7 8 2 8 x 12 - 3 62 114 9 6 6 x 14 + 14 5 4 3 93 x 16 + 12 7 2 3 x 18 + 31 x 20 - 9 x 22 ; 

(4 -x 2 ) V2 h8[x] 2 

a9[x_] = ; 

h7[x] h9[x] 

Reduce [a9[x] > SS < x < 1] 

< x < 1 



The coefficient of sin '(~) in g9 is now constant by this choice of CU). 

f9[x_] = Map[Simplify, Collect [a9 [x] f 8 ' [x] //Together, {si, s2, s3, s4}], 2] 

g9[x__] = Map[Simplify, Collect [a9 [x] g8 ' [x] , ArcSin[_]], 2] 

f9[x] 

r9[x_] = ; 

g9[x] 



7T 2 h9[x] 2 
al0[x_] = ; 

27 648 (4 - x 2 ) 3/2 h8[x] 
Reduce [alO [x] > SS < x < 1] 

< x < 1 
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flO[x_] =Map[Factor, Collect [alO [x] f 9 ' [x] //Together, {si, s2, s3, s4}], 2] 

glO[x_] = alO[x] g9 ' [x] // Together // Factor 

flO[x] 

rlO[x_] = ; 

glO[x] 

We now have f\ a , g i0 , and r i0 all algebraic functions independent of the value of a. 

Reduce [glO[x] >0 SS < x < 1] 

< x < 1 

d[x_] = flO ' [x] glO[x] - flO[x] glO ' [x] // Together; 

dl[x_] =Map[Factor, Collect [Numerator [d [x] ] , {si, s2, s3, s4}], 2] 
d2 [x_] = Denominator [d[x] ] 

Reduce [d2[x] >0 SS < x < 1] 

< x < 1 

dl[x] 

rlO ' [x] == // Simplify 

gl0[x] 2 d2[x] 

True 

Then r[ Q = — — = — , where dj > and g\o > 0; that is, r' l0 > if d[ > 0. 

no * sio 

pi = Coefficient [dl [x] , si] //Factor 
p2 = Coefficient [dl [x] - si pi, s2] //Factor 
p3 = Coefficient [dl [x] - si pi, s3] //Factor 
p4 = Coefficient [dl [x] - si pi, s4] //Factor 

dl [x] == si pi + s2 p2 + s3 p3 + s4 p4 // Simplify 

True 

Further, d\ = i s i Pi is a function rational in the s,-' s. 
Reduce [si pi + s4 p4 > SS < x < 1] 

< x < 1 

Reduce [s2 p2 > SS < x < 1] 

< x < 1 

Reduce [s3 p3 > SS < x < 1] 

< x < 1 

Then d\ > 0, whence r' la > and r m = p 9 S on (0, 1). 

Sign[rlO[0] ] 

1 

So, rio(0 +) > and r )0 S 1 imply r )0 > and hence fi > on (0, 1). Note that this is true for arbitrary a e [0, 1]. This proves Lemma 3.3. 
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TSO: Monotonocity of qr,S;o on (0, 1) 

Following are numerical calculations supporting the arguments of the proof of statement (TSO) in Theorem 2.1. 

{a, b, c} = {0, ARE [0] , ARE ' [0] } 

{0, 1, 0} 

{g9[0], f9[0]} 

{0, 0} 

Sign I® {g8[l] , f8[0]} 

{-1, 1} 

Sign[rl0[0] g8[0] -f8[0]] (* r 9 (0 + ) =r i0 (0 + ) by 1' Hospital, so that p 8 (0 + )<0 *) 

-1 

{g7[0], f7[0]} 

{0, 0} 

Sign/@{g6[0], g6[l/2], g6[3/4], g6[l]} 

{1, 1, -1, -1} 

Sign[f6[0] ] 

1 

Sign[-(r8[0] g6[0] -f6[0])] (* r 7 (0 + ) =r 8 (0 + ) by l'Hospital, so that p 6 (0 + )>0 *) 

1 

p~6 [x_] = - (r7 [x] g6 [x] - f 6 [x] ) ; 
Sign /@ {p~6 [3/4] , p~6[l]} 

{1, -1} 

Sign /@ {g5[0] , g5[l] } 

{1, 1} 

Sign[f5[0] ] 

1 

p~5 [x_] = r6 [x] g5 [x] - f 5 [x] ; 
Sign /@ {p~5 [0] , p~5 [1 / 2] } 

{-1, 1} 

Sign /@ {Limit [g4[x] , x-»0], g4[3/10], g4[2/5], g4[l]} 

{-1, -1, 1, 1} 

Sign /@ {Limit [f 4 [x] , x-»0], f4[2/5], f4[l]} 

{-1, -1, 1} 
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Sign[r5 ' [3/10]] 

1 

p4 [x_] =r5[x] g4 [x] - f 4 [x] //Together; 

ql[x_] = Collect [Numerator [p~4 [x] ] , {ArcSin[_], Sqrt[_]}]; 
q2 [x_] = Denominator [p4 [x] ] ; 

{ql[0], q2[0] }// Simplify 

{0, 0} 

Sign/® {ql ' [0] , q2 ' [0] } 

{0, 1} 

The l'Hospital rule shows that p,(0 +) = - (0 +) = = 0. 

Sign /@ {g3[0] , g3[16/25], g3[l]} 

{0, 1, 1} 

Sign /@ {f3[0] , f3[16/25], f3[l]} 

{0, -1, 1} 

Sign /@ {g2[0] , g2[7/20], g2[l/2], g2[43/50], g2[l]} 

{1, -1, -1, 1, 1} 

Sign /@ {f2[0] , f2[7/20], f2[l/2], f2[43/50], f2[l]} 

{1, 1, -1, -1, 1} 

Sign/@{gl[0], gl[31/50], gl[l]} 

{0, -1, 0} 

Sign /@ {fl [0] , fl[31/50], f 1 [1] /. sub // Simplify} 

{0, 1, 0} 

{g0[0], g0[l]} 

{0, 0} 

{f0[0], fO[l]} 

{0, 0} 

{r2[0] //Simplify, r2[0] / / N} 

3 (32 (-2 + V3~) + ^ 2 ) 
( , 0.0984257) 

1 4tt 2 J 

{r2[l] //Simplify, r2[l] / / N} 

, 135 V3~ - 99 Vl5~ + 40 n -1 

, 0.190467 

L 40 JT J 
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TS1 : Monotonocity of qrr,S;i on (0, 1) 

Following are numerical calculations supporting the arguments of the proof of statement (TS 1 ) in Theorem 2.1. 

{f [1], g[l], f ' [1], g' [1]} / . sub// Simplify 

{0, 0, 0, 0} 

- f"[l] f " ' [1] g" [1] - f " [1] g' " [1] , 

{a, b, c} = 1, , / . sub // Simplify 

g"[l] 3g"[l] 2 1 

9 (-15 + 11 V^) 3 (-135 + 99 -J5~ + 5 (-9 + a/ 5 ") 7r) 
f 1 ' ' 1 

1 40 7T 40 7T 2 ' 

{g9[0], f9[0]} 

{0, 0} 

Sign l@ {g8[l] , f8[0] } 

{-1, 1} 

Sign[rlO[0] g8[0] -f8[0]](* r 9 (0 + ) =r 10 (0 + ) by 1' Hospital, so that p 8 (0 + )<0 *) 

-1 

Sign/@{g7[0], g7[l/5], g7[l]} 

{1, -1, -1} 

Sign[f7[0] ] 

1 

p~7 [x_] = - (r8 [x] g7 [x] - f 7 [x] ) ; 
Sign 19 {pi [0] , pi [1 / 5] , pi [1] } 

{1, 1, -1} 

Sign /@ {g6[0] , g6[l/2], g6[17/20], g6[l]} 

{1, 1, -1, -1} 

Sign[f6[0] ] 

1 

Sign[r7[0] g6[0] -f6[0]] (* p 6 (0 + )>0 *) 

1 

Sign[r7 ' [1/2]] 

-1 

p~6 [x_] = - (r7 [x] g6 [x] - f 6 [x] ) ; 
Sign I® {p~6 [17/20] , P~6[l]} 

{1, -1} 
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Sign /@ {g5[0] , g5[l] } 

{1, 1} 

Sign /@ {f5[0] , £5 [1/2], f5[l]} 

{-1, 1, 1} 

Sign[r6[0] g5[0] -f5[0]] (* p 5 (0 + )>0 *) 

1 



Sign/® {Limit [g4[x] , x-»0], g4[3/4], g4[4/5], g4[l]} 

{-1, -1, 1, 1} 

Sign/® {Limit [f 4 [x] , x-»0], f4[3/4], f4[4/5], f4[l]} 

{1, -1, -1, 1} 

Sign /@ {g3[0] , g3[l/2], g3[l] //Simplify} 

{1, -1, 0} 

Sign /@ {f3[0] , £3 [1/2], f3[l] //Simplify} 

{1, 1, 0} 

Sign/@{g2[0], g2[l/10], g2[l]} 

{-1, 1, 0} 

Sign /@ {f2[0] , f2[l/10], f2[l] //Simplify} 

{-1, -1, 0} 

Sign/® {gl[0] , gl[l]} 

{-1, 0} 

Sign /@ {fl [0] , fl [1] / . sub // Simplify} 

{-1, 0} 

Sign[-(r2[0] gl[0] -fl[0])] (* p 1 (0 + )>0 *) 

1 



{gO[l], fO[l]} 

{0, 0} 



{rO[0] //Simplify, rO[0] // N} 

4 05 - 297 -\[5~ - 6 a/^T (-45 + 19 -y/iT) 7t + 40 jt 



2 

, 0.551626] 

1 4 7T 2 1 

{r4[l] //Simplify, r4[l] / / N} 

1350 (-15 + 11 t/s") - 225 ^75 + a/5~) jt + a/3~ (-1125 + 4297 ^s") " 2 

{ , 1.82004] 

L 2000 7T 3 J 
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RSr: Proof of Lemma 3.4 (Reduction of <fc,s ; a to an algebraic form) 



By default, all functions shall be assumed to be on (0, 1): 

$Assumptions = < x < 1 

< x < 1 

Use the same substitution as given in Appendix TSr: 



sub = {ArcSin [ 1 / 4 ] -» 2 ArcSin [ -\/3/ 8 ] J ; 

The asymptotic mean and variance of Pearson's R: 



Mr[x_] = x ; (* (3.2) in paper *) 
(a 2 ) R [x_] = (l -x 2 ) 2 ; (* (3.1) in paper *) 

The asymptotic mean and variance of Spearman's S: 

6 r x 



Ms [ x _ ] = — ArcSin [—1; (* (3.8) in paper *) 
7T L 2 J 

324 x 2 72 

(a 2 ) [x_]:=l ArcSin[-l + — (II [x] + 2 12 [x] + 2 13 [x] + 4 14 [x] ) ; (* (3.9) in paper *) 

S _2 L q J 7 



ArcSin f , 1 ArcSin f — ^— 1 

L 8-4 X 2 J L 6-2 x 2 J 

II ' [x_] = — — ; 12 ' [x_] 



-x 2 -y^ 



4 - x 



2 



[x 4-x z I r x 4-x z I 
■ I ArcSin I — ' ' I 

2 V 2 " V 8 - 6x2+x4 2 -\/ 1 2-7 x 2 tx 4 

13 ' [x_] = — ; 14 ' [x_] = 



-^4 - x 2 -\J. 



,4 - x 2 

11 [x_?NumberQ] : = NIntegrate [II ' [u] , {u, 0, x}] 

12 [x_?NumberQ] : = NIntegrate [12 ' [u] , {u, 0, x}] 
I3[x_?NumberQ] : = NIntegrate [13 ' [u] , {u, 0, x}] 
14 [x_?NumberQ] : = NIntegrate [14 ' [u] , {u, 0, x}] 
I1[0] =I2[0] =I3[0] =I4[0] =0; 

(a 2 ) s [l] = 0; (* This follows from Lemma 3.1 in the paper. *) 
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Express ARE^ s as a ratio ARE# 5 = -: 



f [x_] := (cr 2 ) s [x] 

36 (1 -x 2 ) 2 



g[x_] 



ARE[x_] = 



7T 2 (4 -x 2 ) 

(° 2 ) S M ^ R '[x] 2 
(a 2 ) R [x] H S ' [x] 2 ' 



f [x] 

ARE [x] == // Simplify 

g[x] 

True 



Begin the reduction phase (cf. (3.11) in the paper): 

f [x_] : = f [x] - b g [x] - c (x - a) g [x] / / Simplify 

gO[x_] = (x-a) 2 g[x] //Simplify 

fO[x] 

rO[x_] = ; 

gO[x] 



al [x_] = yji-x 2 ; 

fl[x_] = Map[Simplify, Collect [al [x] fO ' [x] //Together, {ArcSin [_] } ] , 2] 
gl[x_] = al[x] gO ' [x] // Simplify 
fl[x] 



rl [x_] 



gl[x] 



The S/ below are terms that will be convenient to collect in the functions defined hereafter: 



si = -\^9 - 4 x 2 yjl6 - 12 x 2 + x 4 yjl2 - 8 x 2 + : 

s2 = ^9 - 4 x 2 -^12 - 8 x 2 + x 4 ; 

s3 = yj 9 - 4 x 2 "^16 - 12 x 2 + x 4 ; 

s4 = a/i2 - 8 x 2 + x 4 "\/l6 - 12 x 2 + x 4 ; 



(tl = A^ 



True 



[ t2 -4 



108 - 120 x 2 + 41 x 4 - 4 x 6 I == s2 // Simplify 



True 



[ t3 -4 



144 - 172 x 2 + 57 x 4 - 4 x 6 I == s3 // Simplify 



|t4 = a/± 



True 



sub2 = {tl -» si, t2 -» s2, t3 -» s3, t4 -» s4} ; 
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a2[x_] = (4-x 2 ) 5/2 ; 
f2[x_] = Map[Simplify, 

Collect [Map [Simplify, Collect [a2 [x] fl ' [x] // Simplify // Together / . sub2, {b, c, Sqrt[_]}], 2] /. sub2, 
{b, c, si, s2, s3, s4}] , 2] 

g2[x_] = Map[Simplify, Collect [a2 [x] gl ' [x] , a], 2] 



f2[x] 

r2[x_] = ; 

g2[x] 

At this point, f 2 , g 2 , and hence r 2 are all algebraic functions. 



h3[x_] = x (41 - 20 x 2 + 3 x 4 ) ; 
1 

a3[x_] = ; 

h3[x] 

Reduce [a3[x] > SS < x < 1] 

< x < 1 



With this choice of a 3 , the coefficient of b in / 3 is constant. 



f3[x_] = Map[Simplify, Collect [a3 [x] f2 ' [x] //Together, {b, c, si, s2, s3, s4}], 2] 



g3[x_] = Map[Simplify, Collect [a3 [x] g2 ' [x] , a], 2] 



f3[x] 

r3 [x_] = / . sub; 

g3[x] 



h4 [x_] = 7052 + 30 147 x 2 - 35 490 x 4 + 13 432 x 6 - 2370 x 8 + 189 x 10 ; 

h3[x] 2 

a4[x_] = ; 

h4 [x] 

Reduce [a4 [x] > SS < x < 1] 

< x < 1 



The coefficient of c in f'4 is constant with this choice of 04. 



f4[x_] = Map[Simplify, Collect [a4 [x] f3 ' [x] //Together, {c, si, s2, s3, s4}], 2] 



g4[x__] = Map[Simplify, Collect [a4 [x] g3 ' [x] , a], 2] 



f4[x] 

r4[x_] = ; 

g4[x] 



7T 2 h4 [x] 2 

a5[x_] = ; 

1728 xh3[x] 

Reduce [a5[x] > SS < x < 1] 

< x < 1 



f5[x_] = Map[Simplify, Collect [a5 [x] f4 ' [x] //Together, {si, s2, s3, s4}], 2] 



g5 [x_] = a5 [x] g4 ' [x] / / Factor 



f5[x] 

r5[x_] = ; 

g5[x] 

Now fs, gs, and r$ are all independent of any choice of a. 
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Monotonicity Properties of the ARE: Appendix RS 



Reduce [g5[x] > && < x < 1] 

< x < 1 

d[x_] = f5 ' [x] g5[x] - f5[x] g5 ' [x] // Together; 

dl[x_] =Map[Factor, Collect [Numerator [d [x] ] , {si, s2, s3, s4}], 2] 
d2 [x_] = Denominator [d[x] ] 

Reduce [d2[x] > SS < x < 1] 

< x < 1 

dl[x] 

r5 ' [x] == // Simplify 

d2[x] g5[x] 2 

True 

Hence r' s = — - — "' 5 = -^t, where di > and gs > on (0, 1). 

pi = Coefficient [dl [x] , si] //Factor 

p2 = Coefficient [dl [x] - si pi, s2] //Factor 

p3 = Coefficient [dl [x] - si pi, s3] //Factor 

p4 = Coefficient [dl [x] - si pi, s4] //Factor 

dl [x] == pi si + p2 s2 + p3 s3 + p4 s4 // Simplify 

True 

Further, d\ = 2f =1 i, p, is a rational function in the si s. 
Reduce [pi sl>OSSO<x<l] 

< x < 1 

Reduce [p2 s2>OSSO<x<l] 

< x < 1 

Reduce [p3 s3>OSSO<x<l] 

< x < 1 

Reduce [p4 s4>OSSO<x<l] 

< x < 1 

Then d\ = 2? = i p, > on (0, 1), whence r 5 = p 4 S 1 on (0, 1). 
Sign[r5[0] ] 

1 

As rs(0 +) > 0, one has rs > and hence fs > 0. Note that this fact is independent of the choice of a e [0, 1]. Thus, Lemma 3.4 is proved. 
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Raymond Molzon and losif Pinelis 



RSO: Monotonicity of <7/?,s ; o on (0, 1) 

Following are numerical calculations supporting the arguments of the proof of statement (RSO) in Theorem 2.1. 
{a, b, c} = {0, ARE [0] , ARE ' [0] } 

{o, — , 0} 

L 9 J 

{g4[0], f4[0]} 

{0, 0} 



Sign /@ {g3[x] / . x-> 0, g3[16/25] , g3[l] } 

{-1, 1, 1} 

Sign I® {f3[x] / . x-» 0, f3[16/25] , f3[l] } 

{-1, -1, 1} 

Sign l@ {g2[0] , g2[33/100], g2[l/2], g2[43/50], g2[l]} 

{1, -1, -1, 1, 1} 

Sign /@ {f2 [0] , f2[33/100], f2[l/2], f2[43/50], f2[l]} 

{1, 1, -1, -1, 1} 

Sign/® {gl[0] , gl[3/5], gl[l]} 

{0, -1, 0} 

Sign /@ {fl[0] , fl[3/5], f 1 [1] /. sub // Simplify} 

{0, 1, 0} 

{g0[0], g0[l]} 

{0, 0} 

{f0[0], f0[l] /. sub // Simplify} 

{0, 0} 

{r2[0] //Simplify, r2[0] / / N} 

, 19 8 7 JT 2 , 
+ + , 0.204559 

L 3 A/3~ 36 

{r2[l] //Simplify, r2[l] / / N} 



Monotonicity Properties of the ARE: Appendix RS 



RS1: Monotonicity of qr/?s i on (0, 1) 



Following are numerical calculations supporting the arguments of the proof of statement (RS 1 ) in Theorem 2.1. 



{f [1], g[l], f ' [1], g' [1]} /. sub// Simplify 

{0, 0, 0, 0} 

- f"[l] f " ' [1] g" [1] - f " [1] g' " [1] , 

{a, b, c} = 1, , / . sub // Simplify 

g"[l] 3g"[l] 2 1 

, 9 33 4 , 

1, — + —— , 3- — j 
1 4 4 A /^ 



Sign /@ {g4[0] , g4[4/5], g4[l]} 

{-1, 1, 1} 

Sign/® {f4[0] , f4[4/5], f4[l]} 

{-1, -1, 1} 



Sign I® {Limit [g3[x] , x-»0], g3[l/2], g3[l]} 

{1, -1, 0} 

Sign /@ {Limit [f 3 [x] , x-»0], f3[l/2], f3 [1] // Simplify} 

{1, 1, 0} 

Sign/@{g2[0], g2[l/10], g2[l]} 

{-1, 1, 0} 

Sign/® {f2[0] , f2[l/10], f2[l] //Simplify} 

{-1, -1, 0} 



Sign /@ {gl[0] , gl[l] } 

{-1, 0} 

Sign I® {fl [0] , fl [1] / . sub // Simplify} 

{-1, 0} 

Sign[-(r2[0] gl [0] - f 1 [0] ) ] (* p 1 (0 + )>0 *) 

1 



{gO[l], fO[l]} 

{0, 0} 



{r0[0] //Simplify, r0[0] // N} 

( (l5 -7 a/5 - ) + — ' 0.8 6825 6) 

L 20 * '9 J 

{r4[l] //Simplify, N[r4[l], 8]} (* r (1-) =r 4 (1-) by l'Hospital *) 

, 7 987 



-, 2.6639982| 
4 100 
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